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Abstract 

A Brownian pump in nonlinear diffusive media is investigated in the presence of an unbiased 
external force. The pumping system is embedded in a finite region and bounded by two particle 
reservoirs. In the adiabatic limit, we obtain the analytical expressions of the current and the 
pumping capacity as a function of temperature for normal diffusion, subdiffusion and super diffusion. 
It is found that important anomalies are detected in comparison to the normal diffusion case. 
The super diffusive regime, compared with the normal one, exhibits an opposite current for low 
temperatures. In subdiffusive regime, the current may become forbidden for low temperatures and 
negative for high temperatures. 
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I. INTRODUCTION 



Pumping is an active nonequilibrium transport process in which fluctuations play a very 



important role. Some pumps have already been investigated in the literature 



l|-l6|. Kosztin 



and Schulten [l| studied the fluctuation-driven molecular transport through an asymmetric 
potential pump and three transport mechanisms: driven by potential gradient, by an external 
periodic force and by nonequilibrium fluctuations Moskalets and Buttiker [3] studied heat 
fluxes in adiabatic quantum pumps and developed a approach to the kinetics for an arbitrary 
relation of pump frequency and temperature. Arrachea and coworkers 4] investigated the 
heat transport and the power developed the local driving flelds on an quantum system 
coupled macroscopic reservoirs. They identifled two generic interesting mechanisms: directed 
heat transport between reservoirs induced by the ac potentials and at slow driving, two 
oscillating out of phase forces perform work against each other. Nonadiabatic electron heat 
pump was investigated by Rey and coworkers |5|. They presented a mechanism for extracting 
heat metallic conductors based on the energy-selective transmission of electrons through a 
spatially asymmetric resonant structure subject to ac driving. 

Recently, Sancho and Gomez-marin 6| presented a model for a Brownian pump powered 
by a flashing ratchet mechanism. The pumping device was embedded in a flnite region and 
bounded by particle reservoirs. Their emphasis is on flnding what concentration gradient 
the pump can maintain. All of the above examples have been formulated within a standard 
Brownian framework, for which diffusion properties are normal, that is, with the mean 
quadratic displacement growing linearly with time t. However, it is well know that there 
are media signaled by a growth of the squared dispersion with yU 7^ 1. An important class 
is given by the " porous medium " equation jr iol, which, in the one- dimensional problem 
and in the absence of external forces, can cast in the form 

dtp = Dfid,{p^-'d,p), (1) 

where p is the density of the diffusing substance, x is a dimensionless coordinate representing 

a bond length, angle, or any other chemical or physical state variable, t is the dimensionless 

time, and D is diffusive constant. The mean quadratic deviation follows the law {x{tY) oc 
2 

ti+M. The transport is subdiffusive for yU > 1, normal for p = 1 and superdiffusive for 
< /i < 1. Many physical systems are well-described by this class of processes: percolation 
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of 



^ases in porous media [ll|, dispersion of biological populations 12 1, grain segregation 
13 1 , fluxes in plasma jl^ and nonextensive statistics jl5|. 



The present work is extend to the study of the Brownian pump to the case of the 
anomalous diffusion. We emphasize on finding how particles can be pumped through a 
cell membrane from a particle reservoir at low concentration to one at the same or higher 
concentration in a nonlinear diffusive media. 



II. GENERAL ANALYSIS 



We consider a Brownian pump which is located in a finite region and transports particles 
across the barrier against the concentration gradients by using an unbiased external force. 
Here, we model the device as a pumping system in which overdamped Brownian particles 
moving in a asymmetric finite potential in the presence of an unbiased external force. The 
potential is embedded in a finite region [0,L], 




FIG. 1: Scheme of a Brownian pump: A spatially asymmetric potential ?7o(x) (defined in Eq. (2)) 
is embedded in a finite region and bounded by two particle reservoirs of concentrations po and pi. 
L and A is the length and the asymmetric parameter of the potential, respectively. The particles 
are driven by an unbiased external force F{t), defined by Eq. (3). 



Oji^, \L<x<L, 

where Q is amplitude of the potential and A is its asymmetric parameter. F{t) is an unbiased 
external force and satisfies 

-Fo, riT + ^T <t<{n + l)r. 



where r is the period of the unbiased force and Fq is its magnitude. 

The correlated anomalous diffusion can be described through the following nonlinear 
Fokker-Planck equation js-lOj: 



^p(x,t) = -^[U'{x,t)p{x,t)]+ D^p>^ix,t) 



dx 



(4) 



U{x,t) = Uoix) - F{t)x, (5) 

where (x,t) is a dimensionless 1+1 space time and D = ks is the Boltzmann constant, 
T is the temperature and rj is friction coefficient. The prime stands for the derivative with 
respect to the space variable x. j{x,t) is the probability current and p{x,t) is the particle 
concentration. This nonlinear equation yields anomalous diffusion when /x 7^ 1: subdiffusion 
for /i > 1 and superdiffusion for < yU < 1. 



A. Normal diffusion 



For the normal diffusive case p = 1, Eq. (4) recovers the ordinary Fokker-Planck equation 
and j{x,t) satisfies 

d 

j[x,t) = -U'{x,t)p{x,t) -D—p{x,t). (6) 

The density p{x, t) follows a first order non homogenous linear differential equation, whose 
formal solution is 

p(x,t) = exp[- / ^ ^^^^ dz]{co - 77 / dzexp[f ^ ^^^^ dy]}. (7) 
Jo -U D Jo Jo D 

If F{t) changes very slowly with respect to t, namely, its period is longer than any 
other time scale of the system, there exists a quasistatic state. In the steady state, the 
concentration is iust a function of space thus the fiux becomes a constant ?'. Though, 

n n n 

unlike typical Brownian motors [6|, |16|, 117[ , the boundary conditions are not periodic nor the 
normalized condition is imposed, the unknown constant cq and j can be found by imposing 
the left reservoir concentration po = p(0) and the right concentration pi = p{L) as fixed 
boundary conditions. We can find that Cq = po and 

O{po-piexp|-af]} D F„L 
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where 



Tfp^ r ^ Q-FqXL ^ ,1 , D{X-l)L ^ FoL Q-FoXL 

(9) 

The average current is 

J=- rmt))dt = hj{Fo)+j{-Fo)]. (10) 

T Jo I 

For studying the pumping capacity, we consider the situation in which J tends to zero 
which corresponds the case in which the pump in maintaining the maximum concentration 
difference between the two reservoirs across the barrier with no net leak;ing of particle. This 
situation is analogous the stalling force in Brownian motors. From Eqs. (8-10), we can 
obtain 

_ J(Fo) + /(-Fq) 
Po e^I{Fo) + e— ^/(-Fo) 

B. Anomalous diffusion 

When fi in Eq. (4) is not equal to 1, we have 

jix,t) = -U'{x,t)p{x,t)-D-^p^ix,t). (12) 

In this case, the explicit expression for p{x, t) cannot be extracted from above equation 
by using the same method as the case of /i = 1. Here, we use the method presented by Zhao 
and coworkers 10| to derive the expressions for j and p. Let us assume 

pm{x, t) = G{X, t)p^'{x, t), jm{x, t) = G{x, t)j{x, t) , (13) 



where G{x,t) is a factor function. We assume that there exists the proper pm{x,t) and 
jm{x,t) satisfying diffusion law. 



d 

jmix,t) = -D—pmix,t), (14) 



From Eqs. (13) and (14), we have 



j{x,t) = -Dp^{x,t)-^[\nG{x,t)]-D-^p^{x,t) (15) 



Comparing Eq. (12) with Eq. (15), we have 



G(a:,t)=exp[^^], (16) 



where 

Mx,t) = I U\x,t)ff-\x,t)dx. (17) 

It is obvious that 0^ represents an effective p-dependent potential resuhing from the 
nonhnear diffusive media. For a quasistatic state, we can get the static solution of Eq. (4) 

p,{x) = [^{x)]f~\ (18) 

where [/]+ =max{/, 0} and 

^{x)=p^,-'-^U{x). (19) 
Prom Eqs. (17-19), we can obtain 

0^(a;,t) = -^ln[vl>(a;)]. (20) 
1 /X 



So we can rewrite the Eq. (15) as 



0^^OM) d 0/,(a;,t), 



j{x,t) = -i5exp[-:^:^]-[exp(^:^)p-(a;,t)]. (21) 



Integrating over x from to x, we have 



So we can get the expression for p^{x, t) 

^ - exp[^l • ( 



The constant j can be found by imposing the left reservoir concentration po = p(0) and the 
right concentration pi = p{L) as fixed boundary conditions, 

J(^o) = ^^K-M(Fo)p5'], (24) 



where 
and 

rL 



M{Fo) = [p^o-' + ^FoL]T^., (25) 



/(Fo) = / [^(a;)]T^(ia; 

JO 



^ [(pr^ - AAL) - po] + ^ [(C + SL) - (C + SAL) ] , (26) 
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where A = (-->g;-/-Aa _ g ^ ^"-'f-^^^-m and C = - 

It must be pointed out that when ^(x) is not always positive there exists a cut off of 
probabihty. A cutoff condition (TsaUis cutoff) yields regions with null probability. This is 
because a cutoff of the stationary solution Eq. (18) restricts the attainable space. Since the 
probability of particles visiting the special regions is null, the particles can not pass across 
pump, then there is no current. In order to describe this cut off a function 9{Fq) is defined 

by 



^(i^o) = 

Similarly, the average current is 



1, ^(x) > 0, for all values of x: 

V y ' (27) 

0, otherwise. 



J^\[3{F,)+j{-F,)]. (28) 
Prom Eqs. (24-28), we can obtain the maximum concentration ratio for J = and — 

pI d{Fo)I{-Fo) + d{-Fo)I{Fo) 



^(Fo)/(-Fo)M(Fo) + e(-Fo)/(Fo)M(-Fo) ' 



(29) 



III. RESULTS AND DISCUSSIONS 



Our study focus on the current and the maximum concentration ratio at J = for normal 
diffusion, subdiffusion and superdiffusion. For simplicity, we take /cb = 1, 77 = 1 and L — 1 
throughout the study. 



A. Normal diffusion = 1 

In Fig. 2, we present current as a function of temperature T for normal diffusion [p = 1.0). 
For lower values of pi the current is larger. When pi is more than pq, the current is negative 
for high temperatures. When T — > 0, the particles cannot pass over the potential barrier 
and the current tends to zero. When T — > 00, the ratchet effect disappears, the transport 
is dominated by the concentration difference and the particles move to the left. Therefore, 
there exists an optimized value of T at which the current takes its positive maximum value. 
It is note that the similar behavior is also presented in Ref. 6. 
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FIG. 3: Concentration ratio pi/po as a function of temperature T for J = at Q = 1, A = 0.9, 
Fq = 0.5, and /x = 1. 
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Figure 3 shows the ratio pi /po as a function of temperature T for normal diffusion (// = 1) . 
When T — > 0, no particle can pass over the barrier, thus pi/po ^ oo. As temperature T 
is increased, the ratchet effect reduces and the pumping capacity decreases. Surprisingly, 
the temperature corresponding to the maximum current is not the same the temperature at 
which the concentration ratio for zero current is maximum. This cause for this is that zero 
current induces the maximum concentration ratio, not the minimum one. 



B. Subdiffusion /x > 1 




Figure 4 shows the current J versus temperature T for subdiffusion (p > 1) at pi = 
Po — 1-0. The curve for normal diffusion is observed to be bell shaped. However, as fi 
is increased, the curve becomes to be not smooth. There exist two values of temperature 
at which the curve has inflexions: the lowest temperature to obtain a positive current, the 
optimized temperature to obtain a maximum current. For low temperatures, there lies a 
finite temperature only above which the particle can pass over the barrier, otherwise the 
particles will be confined in both reservoirs. For high temperatures, the current is negative. 
However, when T — > oo, the current will approach to zero. 

9 
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FIG. 5: Currents j(i^o)) —j{—Fo) and J versus temperature T for /x = 2.0 at Q = 1, A = 0.9, 
Fo = 0.5, po = 1.0, and pi = 1.0. The temperatures for the two inflexions are = 0.225, 
rf>^ = 0.725. 



In order to illustrate the transport behavior for subdiffusion, the currents j(-Po)) ~j(~-Po) 
and J as a function of T ior fj, — 2.0 are shown in Fig. 5. Prom Eq. (19) and (27), we can 
find that there exist two values of temperature 

r^^^ = ^(Q±FoA), (30) 
A* 

only above which ^(x) is always positive, otherwise ^(x) may be negative. Therefore, 
when T < T^Fq^ state space becomes disconnected and crossings become forbidden. 
The local current j tends to zero. This is because a cut off of the stationary solution of 
Eq. (18) restricts the attainable space. From the figure, we can see that both j(-Fo) and 
j(— Fq) are zero and J = for T < T|?; j(Fo) is positive, j(— Fq) is zero and J = |j(-Fo 
for Tg <T < T%; j(Fo) is positive,j(-i^o) is negative and J = l\j{Fo) + j{-Fo)] for 
T > T^Fo ■ The first inflexion is at Tg and the second one is at T^Fq ■ 

Figure 6 shows the concentration ratio pi/po as a function of T for subdiffusion {/i — 1.5, 
2.0, and 4.0). The ratio pi/po decreases with increasing T. The curve is not smooth and 
there exists a inflexion at T = Tg. For low temperatures, the ratio pi/po decreases with 
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FIG. 6: Concentration ratio pi/po as a function of temperature T for J = at Q = 1, A = 0.9, 
Fo = 0.5, po = 1-0, and /x = 1.5, 2.0, 4.0. 



increasing /x, while it increases with increasing /j, for high temperatures. 
C. Superdiffusion < /x < 1 

Figure 7 shows current J as a function of temperature T for superdiffusion (0 < // < 1) 
at pi — Po- It is found that the curve is not smooth and there is a inflexion at which the 
negative current takes its maximum value. For low temperatures, the particles move to the 
left and the negative current increases with decreasing fi. As temperature T is increased, 
the current becomes to be positive and the positive current increases with decreasing p,. 

In order to explain the transport behavior for superdiffusion, the currents j{FQ), j{—Fo) 
and J versus temperature T ior p — 0.1 are shown in Fig. 8. We can find from Eqs. (19) 
and (27) that there is a critical value of temperature 

Tg = -Fof^, (31) 
P' 

below which the local current j(-Fo) is zero. When T is less than Tg, j{Fo) is zero and 
j{—Fo) is negative, so the current J = |j(— Fq) is negative. When T is more than Tg, j{Fo) 
is more than — _;(— Fq), so the current J is positive. As temperature T is increased, j{Fo) 
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FIG. 7: Current J as a function of temperature T for different values of /Lt at Q = 1, A = 0.9, 
Fo = 0.5, po = 1-0, and pi = 1.0. 




FIG. 8: Currents 3{Fq), j{—Fo) and J versus temperature T for /x = 0.1 at Q = 1, A = 0.9, 
Fo = 0.5, po = 1-0) and pi = 1.0. The temperature at the inflexion is TS = 4.5. 



12 



tends to — j(— Fq), then the current J goes to zero. 




FIG. 9: Concentration ratio pi/po as a function of temperature T for J = at Q = 1, A = 0.9, 
Fo = 0.5, po = 1.0, and p = 0.5, 0.6, 0.7, 0.8, 0.9. 

Figure 9 shows the concentration ratio pi/po as a function of temperature T for superdif- 
fusion {n — 0.5, 0.6, 0.7, 0.8 and 0.9). For low temperatures, T is less than and j(-Fo) 
goes to zero, so the particles move to the left and pi/po < 1. As is decreased, the lowest 
temperature for pumping the particles to the right increases. Therefore, the pumping system 
in superdiffusive regime need a higher temperature than that in normal diffusive regime. 

IV. CONCLUDING REMARKS 

In this study, we investigate a Brownian pump in nonlinear diffusive media with an 
unbiased external force. The pump is embedded in a finite region and bounded by two 
particle reservoirs. The analytical results arc obtained in the adiabatic limit. In normal 
diffusive regime, current is a peaked function of temperature and concentration ratio for zero 
current decreases with increasing temperature. In subdiffusive regime, current is forbidden 
for low temperatures and negative for high temperatures. The concentration ratio for zero 
current decreases with increasing temperature, also. In superdiffusive regime, the transport 
compared with normal diffusion exhibits an opposite direction for low temperatures. There 
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exists a finite temperature at which the concentration ratio takes its maximum value. For 
anomalous diffusion, there exist inflexions in the curves for J vsT and pi/po- This is because 
a cut off of the stationary solution Eq. (18) restricts the attainable space. 
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